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Abstract 

The a 2 s correction to the Ellis-Jaffe sum rule for the structure function g± of polar- 
ized deep inelastic lepton-nucleon scattering is calculated. 
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The results of the EMC Collaboration at CERN and E130 Collaboration at SLAC 
[0] for the Ellis- Jaffe || sum rule Jq dxg\(x, Q 2 ) attracted a lot of attention to this sum 
rule; see HI- [Oil and references therein. Recent data of the SMC Collaboration at CERN 



HJ on polarized scattering of muons off deuterium and of the E142 Collaboration at SLAC 
14] on polarized scattering of electrons off helium 3 He allowed the determination of the 



analogous sum rule Jq 1 dxg™(x, Q 2 ) for a neutron. This in turn allowed us to find a difference 
Jq 1 dx[gi(x, Q 2 ) — <7i(x, Q 2 )} which is the Bjorken sum rule ][HJ. At present, the Bjorken sum 



rule is calculated within QCD with quite high accuracy. The a s correction [16|, a s correction 



17] and a z s correction |1| are calculated in the leading twist approximation. The higher 



twist corrections are also calculated [|T{J. For the Ellis- Jaffe sum rule only the a s correction 
was calculated in the leading twist [p0[] . The power corrections were calculated in [[2l|j . 

In the present paper we obtain the a 2 correction to the Ellis- Jaffe sum rule in the leading 
twist massless quark approximation. All calculations are performed in dimensional regular- 
ization 22l. Renormalizations are done within the MS'-scheme II23U , the standard modifica- 



tion of the Minimal Subtraction scheme fl2~3 



Polarized deep inelastic electron-nucleon scattering is described by the hadronic tensor 



W, u (p,q) = — 
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of which we will consider the structure function g\. Here = ipj^Eip 



the electromagnetic quark current and E 
magnetic charge matrix, x 



Zi^iln^i is 

diag(2/3, —1/3, —1/3, ...) is the quark electro- 

q 2 . The nucleon state 
^2j9°(2vr) 3 (5( 3 )(p - p'). s a is the 



is the Bjorken variable, Q 2 



| p, s > is covariantly normalized as < p, s | p , s > 
polarization vector of the nucleon: s a = U(p, s) r y fT '~f5U(p, s), where U is the nucleon spinor, 
U{p,s)U(p, s) = 2M. 

The moments of the deep inelastic structure functions are expressed |25| via quantities 
of the Wilson operator product expansion (OPE) of the corresponding currents. We need 
the OPE of two electromagnetic currents. The strict method of the OPE ensures p6|,|27|| 
that the OPE of two gauge-invariant currents can contain only gauge-invariant operators 
with their renormalization basis. Thus we have only contributions from the non-singlet and 
singlet axial currents in the OPE of electromagnetic currents in the leading twist for the 
considered structure: 

'dz e^T{J^z).M0)} Q2 =°° 
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J2 C a (\og(^-), a s (fi 2 ))J^(0) + C s (\og(^-), a s (^))4(0) + higher twists 



(2) 

where the non-singlet contribution can be rewritten as 

^CV ff 5 ' tt (0) = C ns Y,Tr{E 2 t a )Jl> a {0), 

a a 

to introduce, as usual, the unique non-singlet coefficient function C ns not depending on the 
number a. 

Here J^ a (x) = ip(x)^( cr r )^t a ip(x) is the non-singlet axial current, where t a is a generator 
of a flavour group, Tr{t a t b ) = ^S ab , and J„(x) = J2ili' l Pi{. x )'lal5' l Pi(x) is the singlet axial 
current. The known twist-two and spin-one axial gluon current K a = 4e aUlU2U . i (Al i d U2 Al 3 + 
\gf ahc A a Vx A h V2 Al^ has also the necessary quantum numbers, but it cannot contribute to 
the above operator product expansion because it is not gauge invariant. Here and further 
on (before presenting the final results), we use the most practical definition for the strong 
coupling constant from the calculational point of view 
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The Ellis- Jaffe sum rule is expressed as 



£ dxg? n \x,Q 2 ) = C ns (l,a s (Q 2 )) {±^\9a 

( ra s (Q 2 ) >y s (n')\ 1 



(3) 



where some comments are in order. Here p(n) denotes a target: proton (or neutron). The 
plus (minus) before \g A \ corresponds to the proton (neutron) target. The proton matrix 
elements of the axial currents are defined as follows: 



g A \sa = 2 < p, s | J^' 3 \p,s>= (Au - Ad)s a , 



where g A / 9v = —1-2573 ± 0.0028 |28| is the constant of the neutron beta-decay; 

a 8 s a = 2V3 <p,s\ J^' 8 I p, s >= (Am + Ad - 2As)s CT , 

S(/i 2 )s CT =<p,s\4\p,s >= (Au + Ad + As)s a , 
and we use the standard notation 

Aq(fi 2 )s a =< p, s I qialtf I p, s >, q = u,d, s. 

We omitted the contributions of the nucleon matrix elements for quarks heavier than the 
s-quark but it is straightforward to include them. We should stress here that g A and ctg do 



2 



not depend on the renormalization point /z 2 since the corresponding non-singlet currents J^' 3 
and J^' 8 are conserved in the massless limit, and hence their renormalization constants are 
equal to one. On the contrary, the singlet axial current has a non-trivial renormalization 
constant. Hence the quantity £(/z 2 ) does depend on the renormalization point (that is why 
it is not a physical quantity). The coefficient functions C ns (l, a s {Q 2 )) and C s (l,a s (Q 2 )) are 
normalized in the standard way to the unity at the tree level. The renormalization group 

2 

technique was applied to the coefficient functions to kill logarithms log(^). 

The singlet axial current has the non-zero anomalous dimension '~y s (a s ) due to the axial 
anomaly []29, |30|. So one can say that the axial anomaly contributes to this sum rule through 



the renormalization group exponent of the singlet current contribution in eq.([|). 

Let us define the functions 7 s (a s ) and (3(a s ) in the renormalization group exponent of 
eq.(|3|). The renormalization group QCD /5-function is calculated |32| in the MS-scheme 
at the 3-loop level : 



l2 ~d^ = _ ^° a * ^ ia * ^ s 

2 \ „ / 38 \ , Z2857 5033 325 
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To define the anomalous dimension of the singlet axial current we should first define the 
singlet axial current itself within dimensional regularization. To define the singlet axial 



current we will follow the lines of ref. |]33fl , where the 't Hooft -Veltman definition [22] of 
the 75-matrix is elaborated for the multiloop case. The singlet current J\ is renormalized 
multiplicatively and is expressed via the bare one [J^]b as 

4 = ZlZ s MS [Jl\ B . (5) 



Here Z S MS is the MS renormalization constant which contains only poles in the regular- 
ization parameter e, the dimension of the space-time being D = 4 — 2e. The extra finite 



renormalization constant Zf is introduced to keep the exact 1-loop Adler-Bardeen form [33] 
for the operator anomaly equation within dimensional regularization in all orders in a s : 



d,4 = a s ^(GG), (6) 



where all quantities are renormalized ones. GG = e ^A/?G" G". and G* = d^A^ — d u A a ^ + 
gf abc A b ^A c v is the gluonic field strength tensor. 

In fact the full physical quantity, the Ellis- Jaffe sum rule, does not depend on the choice 
of the normalization constant Z§. But to be definite we adopt the normalization of [53] to 



keep in the MS-scheme the singlet axial current satisfying eq.(^|). 

The anomalous dimension of the singlet axial current is zero at the 1-loop level and 
starts from the 2-loop level. To have the next-to-leading approximation we need two non- 
zero terms, i.e. the 3-loop approximation. The 3-loop approximation for the anomalous 
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dimension of the singlet axial current was calculated in [[!3j and confirmed in The 
result in the adopted normalization reads 



7 s (a s ) = /i ; 
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3 / J 3 

3 are the Casimir operators of the fundamental and adjoint rep- 
resentation of the colour group SU(3). The a 2 term agrees with the calculation |20| after 
multiplication of our result by the factor (—2) due to different normalizations. 

The non-singlet coefficient function C ns (l, a s {Q 2 )) was calculated in the approximation 



in 



18| where the Bjorken sum rule was calculated in this approximation. We want to obtain 



the a 2 correction to the singlet contribution to the Ellis- Jaffe sum rule ((^): 

C s (l,a s (Q 2 ))exp 



= C s (l,a s (Q 2 )) l-a s (Q 
where we introduced the notation 
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for the renormalization group-invariant (i.e. /^-independent) nucleon matrix element of the 
singlet axial current. 

Beside the 3-loop (the order a^) approximation of the anomalous dimension, we need also 
the 2-loop (the order a 2 ) approximation for the singlet coefficient function C 8 (l, a s ) which has 



already been calculated in | 36fl . Here we present the calculation of C s with another method 



to confirm the validity of the result obtained in [3H]. We use the "method of projectors' 
To project out the coefficient function C s from the OPE of eq.(f^), one should sandwich this 
equation between quark states and nullify the quark momentum p. To be more precise, one 
should consider the following Green function 



, I ,Jz e** < | Til>(ph*'YBi/>(j>)Jr(z)J u {0) | >|^ ted . 



w^^(iog(fe)> <» 2 )) < o I m P ) lal ^( P )z s 5 z s MS [4(o))]B I o >i;x 4ated , (io) 

where some remarks are in order. ip(p) is the Fourier transform of the quark field carrying the 
momentum p. Quark legs are amputated. The essence of the method is the nullification 
of the quark momentum p. In the dimensional regularization scheme all massless vacuum 
diagrams are equal to zero. So on the r.h.s. only the tree graphs survive after the nullification 
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of p. In our case the only operator which produces a tree graph is J%. The non-singlet axial 
current does not contribute because of the nullification of the flavour trace: Tr(t a ) = 0. It is 
interesting to note that at p = we have infrared divergences in the diagrams of the l.h.s. of 
eq. (ITHD . But these divergences are cancelled by the ultraviolet poles of the renormalization 
constant Z S MS of the singlet current. 

Thus to calculate C s we need to calculate the diagrams contributing to the l.h.s. of 
eq . (|T0|) . These are the diagrams of the forward scattering of a photon off quarks with 
photon momentum q and zero quark momentum. In comparison with the calculation of the 
non-singlet coefficient function C ns |17|, |18 , we have at the 2- loop level two extra diagrams 



where both electromagnetic vertices are inside a closed quark loop. The analytic calculation 
of the diagrams has been done with the symbolic manipulation program FORM by 



means of the package MINCER p9 ]. This package is based on algorithms of ref. ||40 |. The 
result is 



C s (l, a s {Q 2 )) = l+a s {-3C F ) + a] 
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, - 23C F C A + ( 8C 3 + y ) C F n f 



(11) 



where (3 is the Riemann zeta-function (£3 = 1.202056903...). This result of ours agrees 
with the calculation in if one takes into account the fact that another finite constant Z£ s 
(relevant for the non-singlet axial current, see [|I8| , |33| ) was taken in for the normalization 
of the singlet axial current instead of our Z\. Multiplying our result ( |TTD by the factor 

a 2 s (22C 2 F - ±fC A C F + %C F n f ) 



yns 



1 + a s (-AC f 
1 + a(-4C F ) + a\22C 2 F - ^fC F C A + |C F n/) 



+ 0(ai 



one can reproduce the result of p6| . So we have a strong check of C s . 

In principle our technique allows us to compute also the ct^ correction to the singlet 
coefficient function C s . But we need then also the 4-loop anomalous dimension of the 
singlet axial current in order to take into account this correction self-consistently in the 
next-next-to-leading approximation for the the Ellis- Jaffe sum rule. The calculation of the 
4-loop anomalous dimension is very time-consuming at present, although all the necessary 
techniques are available. 

Collecting together all the relevant results for the coefficient functions and the anomalous 
dimension, we obtain finally the next-to-leading approximation for the Ellis- Jaffe sum rule: 
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^2 for the strong coupling constant. We keep the known (extra for 
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Here we use — 

7T 

the next-to-leading approximation) term |TB] for the non-singlet part, since this part 
determines the Bjorken sum rule. For the singlet part we factorize the Q 2 -dependent fac- 
tors coming from the coefficient function (the first factor in the singlet part) and from the 
renormalization group exponent (the second factor). The leading a s term agrees with pO 
For the case rif = 3 the sum rule reads 
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The Ellis- Jaffe sum rule looks most compact for the choice fi 2 = Q 2 since the renor- 
malization group exponent becomes a unity. But then the Q 2 -dependence jumps inside the 
matrix element of the singlet axial current X: 

>2\\ 2 /„, ^2\\ 31 



f 1 dxgl {n \x,Q 2 
Jo 



a s (Q 2 



3.5833 



a s (Q 2 



20.2153 



a s (Q 2 



7Y 



7Y 



7Y 



, 1 , , 1 

x ±— \g A \ + 



"12 1 



36 



+ 



'a s {Q 2 



1.0959 



'a s {Q 2 



Ti 



7i 



9 



S(Q 2 ). 



(14) 



The conclusion is that the a 2 correction to the Ellis- Jaffe sum rule is quite small in the 



MS'-scheme. 
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